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1.  INTRODUCTION 


Test  range  measurements  for  fired  artillery  shells  are  of  great  importance  to  ballisticians,  members 
of  the  scientific  community  who  study  the  motion  of  artillery  projectiles  and  missiles.  The  accuracy 
of  aiming  data  developed  by  ballisticians  can  only  be  expected  to  be  as  good  as  the  accuracy  of  test 
range  measurements  taken  at  a  proving  ground.  It  follows  that  an  understanding  as  to  how 
measurements  of  range,  deflection,  time  of  flight,  and  height  of  burst  for  a  fired  artillery  shell  are 
obtained  should  be  of  value  to  the  ballistician.  Actually,  such  knowledge  would  benefit  any  weapon 
systems  engineer  who  must  make  decisions  either  directly  or  indirectly  by  depending  on  the  use  of 
aiming  data. 

This  report  will  deal  with  the  actual  engineering  upon  which  triangulation  is  based,  as  well  as  the 
analytical  and  mathematical  procedures  by  which  observer  measurements  are  converted  into 
measurements  of  range,  deflection,  height  of  burst,  and  time  of  flight.  It  is  designed  to  give  an 
individual  an  introduction  to  the  triangulation  technique.  However,  a  copy  of  the  triangulation 
analytical  computer  program  is  included  (Appendix  B),  as  well  as  sample  data  input  (Appendix  C)  and 
corresponding  data  output  (Appendix  D).  This  will  also  allow  one  to  study  triangulation  in  as  much 
detail  as  desired. 


2.  THE  OPERATIONAL  ASPECTS  OF  TRIANGULATION 

Suppose  one  wants  to  obtain  test  range  measurements  for  a  M483A1  projectile  fired  at  a  quadrant 
elevation  of  950  mils  from  a  Ml 98  howitzer  using  a  particular  charge.  Also,  suppose  that  one  wants 
the  projectile  to  burst  at  a  height  of  approximately  200  meters  above  ground. 

The  first  step  is  to  estimate  how  far  the  projectile  will  fly  and  to  determine  the  fuze  setting 
required  to  obtain  the  desired  height  of  burst.  From  a  tabular  firing  table  or  field  computer,  the  test 
engineer  will  determine  approximate  range,  deflection  (drift),  and  fuze-setting-to-burst  values.  Then, 
knowing  the  position  coordinates  of  the  howitzer  and  its  azimuth  of  fire  from  survey,  the  test  engineer 
will  have  an  idea  as  to  where  down  range  the  projectile  can  be  expected  to  burst. 
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Accordingly,  ju  to  four  observers  are  stationed  in  towers  most  opportune  in  location  to  the 
expected  burst  and  impact  areas.  From  survey,  the  location  and  heights  of  the  towers  are  known. 

Each  tower  observer  will  make  use  of  an  instrument  called  a  modified  transit,  which  is  a  transit 
equipped  with  a  monocular  lens  with  crosshairs  that  magnifies  and  has  a  wide  field  of  view.  With  this 
instrument  each  observer  can  measure  azimuth  and  elevation  angular  values  with  respect  to  the  tower 
after  sighting  the  burst  event  through  the  monocular  lens. 

During  the  shoot,  each  observer  will  be  in  radio  contact  with  the  gun.  Thus,  upon  firing,  each 
observer  will  know  when  to  start  his  stopwatch  and  will  know  when  to  expect  burst  or  impact. 
Whenever  video  is  not  used,  it  is  customarily  the  responsibility  of  the  observer  to  obtain  time  of  fli.  vt 
of  the  round  fired  as  well  as  angular  measurements  to  burst  or  impact. 

In  a  test  situation,  practice  rounds  called  spotters  are  usually  shot  first  so  that  it  becomes 
relatively  certain  that  each  observer  will  be  successful  in  obtaining  all  measurements  for  the  actual 
firing  of  the  round. 

It  follows  from  the  engineering  behind  triangulation  that  the  test  range  measurements  obtained  for 
a  fired  projectile  are  as  follows: 

1)  tower  locations  and  heights 

2)  azimuth,  elevation,  and  time  of  flight  measurements  from  each  tower 

3)  location  and  azimuth  of  gun. 

The  test  range  measurements  obtained  above  must  be  converted  into  range,  deflection,  height  of 
burst,  and  average  time  of  flight.  The  next  section  of  this  report  will  show  how  mathematical 
procedures  can  be  used  to  do  this. 
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3.  TRIANGULATION 


Triangulation  is  the  analytical  mathematical  procedure  for  obtaining  ranges,  deflections,  heights  of 
burst,  and  time  of  flights  for  fired  projectiles  from  raw  test  range  data.  For  an  air  burst  round,  one 
first  wants  to  compute  the  range  and  deflection  from  the  gun  point  to  the  perpendicular  projection  of 
the  air  burst  point  to  the  ground  ("impact  point").  The  range  is  the  distance  from  the  weapon  to  the 
perpendicular  projected  point  on  the  ground  and  the  deflection  is  the  angular  and  linear  distance  from 
the  weapon  line-of-fire  to  the  same  point. 

"Impact  points”  are  computed  from  azimuth  readings  taken  from  one,  two,  three,  or  four  towers. 
There  arc  four  modes  of  computing  impacts  -  one  for  each  tower  combination. 

For  one  tower,  the  "impact  point"  is  determined  by  the  intersection  of  the  gun  azimuth  and  tower 
azimuth.  On  the  computer  output,  the  deflection  will  be  dashes.  The  range  is  the  distance  from  the 
point  of  intersection  to  the  weapon  jXisition  (Figure  1). 

For  two  towers,  the  "impact  pemt"  is  determined  by  intersecting  both  tower  azimuths.  Range  is 
defined  as  above.  Deflection  from  linc-of-firc  is  computed  as  a  perpendicular  distance  from  "impact" 
to  the  line-of-fire  and  as  an  angle  between  the  linc-of-fire  and  vector  from  weapon  to  "impact  point" 
(Figure  2). 

For  three  towers,  the  "impact  point"  is  computed  as  the  average  of  three  intersection  points  of  the 
three  azimuth  readings  with  the  restriction  that  the  area  of  the  triangle  determined  by  the  three 
intersection  points  be  less  than  a  specified  test  area,  usually  50  square  meters.  If  the  area  is  greater 
than  this  specified  test  area,  the  azimuths  are  considered  erroneous  and  the  output  consists  of  azimuth 
corrections.  If  the  computed  area  is  less  than  or  equal  to  the  test  area,  the  "impact”  is  taken  to  be  the 
average  of  the  three  intersection  points.  The  range  and  deflection  are  computed  in  the  same  manner 
described  in  the  case  for  two  towers  (Figure  3). 

For  four  towers,  the  towers  are  grouped  three  to  a  set  to  form  various  combinations.  That 
combination  which  forms  the  smallest  triangular  area  will  be  used  to  determine  the  range  and 
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deflection,  provided  its  area  is  less  than  or  equal  to  the  specified  test  area.  "Impact"  is  the  average  of 
the  three  intersections.  Range  and  deflection  arc  computed  as  before  (Figure  4). 

The  geometry  and  trigonometry  of  how  range  and  deflection  are  computed  for  two,  three,  or  four 
tower  readings  are  illustrated  in  Figure  5. 

Computing  altitude  (height  of  burst)  for  an  airburst  round  from  the  elevation  measured  from  a 
tower  is  illustrated  in  Figure  6.  Altitudes  arc  computed  for  each  tower  and  then  averaged.  Times  of 
flight  measured  with  stop  clocks  from  each  tower  are  also  averaged. 

The  triangulation  computer  program  written  for  the  Aberdeen  Proving  Ground  (APG)  by  Dennis 
Flaherty  and  related  information  are  shown  in  Appendices  A,  B,  C,  and  D.  It  should  be  noticed  in  his 
program  that  azimuths  are  measured  clockwise  from  due  south.  Notice  also  that  the  cartesian 
coordinate  system  with  origin  at  the  gun  has  its  x-axis  inclined  35  degrees  from  due  south  and  is  right- 
handed  with  the  positive  z-axis  going  into  the  ground.  The  program  can  easily  be  modified  for  use  at 
Yuma  and  Dugway  Proving  Grounds  where  azimuths  arc  measured  clockwise  from  due  north  and 
where  the  cartesian  coordinate  system,  with  origin  at  the  gun  has  its  x-axis  due  north  and  its  y-axis  due 
cast. 


4.  DISCUSSION 

The  accuracy  of  triangulation  is  very  good  if  we  can  assume  that  at  least  three  tower  azimuth 
readings  arc  good  and  if  we  can  assume  that  the  true  impact  is  within  the  formed  triangle.  Making 
such  assumptions  and  supposing  that  the  formed  triangle  is  a  square  right  triangle  of  50  square  meters, 
the  range  and  deflection  errors  cannot  be  larger  than  7.46  meters. 

For  small  ground  impact  firing  tests,  it  is  possible  for  a  person  to  mark  the  crater  location  formed 
by  the  impacting  round  with  a  stake  and  to  have  all  towers  triangulate  on  the  stake.  If  one  does  this, 
one  should  obtain  near  perfect  accuracy  for  range  and  deflection  values.  However,  this  process  is  very 
laborious  and  time  consuming,  and  is  impractical  to  be  used  for  large  test  Firing  plans.  In  such  cases, 


4 


and  for  air  burst  occasions,  burst  sightings  from  towers  are  relied  upon.  When  so,  it  is  conceivable 
that  all  towers  can  triangulate  perfectly  but  because  of  like  systematic  tower  errors,  be  off  the  true 
mark. 


Triangulation,  which  is  based  on  burst  sightings,  does  have  its  critics.  Mr.  Robert  Lieske,  a 
supervisor  and  colleague  of  mine,  has  argued  that  when  four  tower  readings  are  used  under  such 
circumstances  and  all  four  readings  are  relatively  good,  triangulation  does  not  use  all  the  information  it 
has  access  to.  When  using  many  towers,  he  has  suggested  an  alternative  technique  based  on  a  least 
squares  method  which  minimizes  angular  error  and  which  selectively  uses  all  the  information. 

Mr.  Lieske  believes  that  minimizing  angular  error  is  more  efficient  than  minimizing  the  area  of  a 
triangle.  Hopefully,  this  technique  can  be  explained  in  a  subsequent  report. 


5.  CONCLUSION 


In  conclusion,  it  is  hoped  by  the  author  that  the  objectives  in  writing  this  report  and  mentioned  in 
the  introduction  have  been  achieved.  The  triangulation  computer  program  is  on  the  Vax  8600  and  780 
computers  in  the  Launch  and  Flight  Division  of  BRL. 
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Figure  1.  One  Tower  Reading. 
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Figure  3.  Three  Tower  Readings. 
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Figure  4.  Four  Tower  Readings. 
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LISTING  OF  INPUT  TO  TRIANGULATION  PROGRAM 


All  input  is  punched  into  a  file,  using  the  following  format. 
For  IRX  =  01  or  02 


1  terns 

Forma  t 

IRX  (01  or  02) 

GUN 

XLOC  ( location ) 

DATE 

PRENG  (Project  Engineer) 

ELEV  (Elevation  -  ft.),  ABOVE 

KRD ,  N ,  NSTEP,  A 

12 

14A4 

13A4 

14A4 

11A4 

F6.2,  IX 
314,  F20 

X  V 

1  1 

N  <  4 

2F11 .3 

x  y 

n  n 

X  ,  Y  ,  AZD,  AZM 

5  5 

4F11.3 

For  N  of  Towers: 

AZD ,  AZM,  ALX,  BEX,  ALPHA,  TIME,  SID  2F7.2,  2F7.2,  A8 ,  F6.2,  A4 

If  blank  data  use  -999  for  both  AZD  and  ALX 
If  blank  data  for  time,  leave  TIME  blank  on  input  card 


AZD 

n 

+  999. 


For  IRX  =  03  or  04 

IRX  12 

♦U,  V,  HI,  HCOR  4F11.3 

♦For  each  tower 

IRX 

999999. 

05  Termination 

Note:  Maximum  number  of  tower  readings  199. 

Note:  When  using  IRX  of  03  be  sure  to  use  specified  point  coordinate 

for  gun  coordinates  if  no  gun  coordinates  are  given. 
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DEFINITIONS  OF  INPUT  TO  TRI ANGULATION  PROGRAM 


IRX  An  integer  e  _,ual  to  1,2,3,  or  4  -  it  specifies  the  type  of  run. 

1  -  Compute  range  and  deflection 

2  -  Computes  impact  coordinates  in  APG  grid  system 

3  -  Computes  distance  from  specified  point  (e.g.  tower  or 

weapon)  to  impact  point 

4  -  Computes  altitude 

GUN  56  character  alphanumeric  field  weapon  identification 

XLOC  56  character  alphanumeric  field  weapon  location 

DATE  56  character  alphanumeric  field  date 

PRENG  44  character  alphanumeric  field  project  engineer 

ELEV  Height  of  weapon  (feet) 

ABOVE  24  character  alphanumeric  field  -  specifies  reference  level 
for  weapon  elevation 
KRD  sequence  number  of  first  round 

N  number  of  towers 

NSTEP  Increment  for  round  number 

A  Limiting  area  (used  with  3  or  4  towers) 

X  Y  Position  coordinates  of  i-th  tower  (i  4)  in  APG  grid  system 

i  i 

X  Y  Position  coordinates  of  weapon,  APG  system 

5  5  ' 

A2D ,  AZM  Line  of  fire  measured  from  south,  in  degrees  and  minutes 
AZD,  AZM  (second  set)  Tower  readings  measured  south  in  degrees  and 
minutes 

ALX ,  BEX  Tower  elevation  readings  of  burst  height,  in  degrees  and 
minutes 

ALPHA  8  character  alphanumeric  field  used  to  label  a  tower  reading 

TIME  Observed  time  of  flight  from  a  given  tower,  in  seconds 

SID  4  character  alphanumeric  character  S  to  be  used  to  identify 

?  spotter  round  as  opposed  to  a  numbered  round 
U,  V  Position  coordinates  of  tower  (used  for  IRX  *  03  or  04)  in 

APG  system  (meters) 

HI  Elevation  correction,  in  meters  (used  for  IRX«*04) 

HCOR*  Elevation  correction,  in  minutes  (used  for  IRX»04) 

*As  of  26  June  75,  there  will  no  longer  be  a  requirement  of  HCOR  for 
APG  firings. 
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DESCRIPTION  OF  OUTPUT  TO  TR I ANGULATI ON  PROGRAM 


Taking  the  various  conditions  in  order,  the  output  formats  are 
as  follows: 

I.  IRX  «=  1 

Output  is  stored  until  the  completion  of  the  processing  run.  If  any 
error  rounds  are  encountered,  the  checking  routine  prints  out  the  correc¬ 
tions  as  follows: 

RANGE  FIRING  SUMMARY  ERROR  CORRECTIONS 

Round  No.  AREA  (16,  D16.9)  The  area  has  been  calculated  as  larger  than 

the  standard  value  given  on  input 
Tower  X  Azimuth  X  (  2 1 4 ,  IX,  F7.3) 

Tower  Y  Azimuth  Y 

Tower  Z  Azimuth  Z  To  correct  such  an  error  one  of  the  tower  readings 

is  changed  to  the  azimuth  given  or  the  standard 
The  azimuths  are  given  in  area  is  increased  to  be  greater  than  the  calculated 
integral  degrees,  and  min-  area, 
utes  to  three  decimal  places. 

If  no  errors  are  encountered  during  the  processing  run  the  normal  type  output 
process  is  reached.  A  page  would  resemble  the  following: 

RANGE  FIRING  SUMMARY  COPY  1  OF  3 

GUN  175  MM  GUN  LOCATION  X  RANGE 

DATE  28  MAY  1969  PROJECT  ENGINEER  MR.  SMITH 

ELEVATION  OF  TRUNNION  6.73  FT  ABOVE  MLW 


AZIMUTH  LINE  OF  FIRE 


RDS  1  to  20  AZ 

35D 

3.72  M 

Rd 

Deflection 

Range 

Times  of 

Flight 

No. 

mils  meters 

meters  Ord  of 

De  t 

1  2 

3 

4 

Mean 

16 

F6.1  F6.1 

F8.1  A8 

F6.2  F6.2 

F6.2 

F6.2 

F6.2 

Three  copies  are  printed  for  each  run. 
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IRX 


2 


II . 

The  error  checking  procedure  and  output  are  the  same  as  for  IRX  -  1. 

The  standard  output  is  as  follows. 

RANGE  FIRING  SUMMARY  COPY  1  of  3 

GUN  175  MM  LOCATION  X  -  RANGE 

DATE  28  MAY  1969  PROJECT  ENGINEER  MR.  SMITH 

ELEVATION  OF  TRUNNION  6.73  ft  above  MLW 

AZIMUTH  LINE  OF  FIRE 

Rds  1  to  20  AZ  350D  3.72  M 

GRID  COORDINATES  OF  POINT  OF  IMPACT 

ROUND  X  Y 

16  F17.2  F17.2 

Again  three  copies  are  provided. 

III .  IRX  -  3 

RANGE  FIRING  SUMMARY  COPY  1  of  3 

ROUND  N*  (RANGE,  TANGENT  BHT )  AV  BHT  DIFF 

16  N  (IX,  F8.1,  IX,  F9 . 5 ,  IX,  F5 . 1 )  F5.1  F5.1 

N  represents  the  number  of  towers  or  other  reference  taking  readings 
1  N  4 

ROUND  AVBHT  DIFF  These  are  copies  of  the  respective  columns 

16  F5.1  F5.1  above.  Only  one  copy  of  this  set  is  provided. 

There  are  three  for  the  above. 
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APPENDIX  B: 

LISTING  OF  TRIANGULATION  COMPUTER  PROGRAM 
WRITTEN  FOR  THE  ABERDEEN  PROVING  GROUND,  MARYLAND, 
BY  DENNIS  FLAHERTY 
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c  PROGRAM  =  52.00  0001 

C  IMPACT  LOCATION  FROM  TOWER  OBSERVATIONS  1,  2,  3.  OR  4  TOWERS  0002 

C  FORTRAN  IV  0003 

C  APRIL  1966  0004 

C  PROGRAM  REVISED  MARCH  1971  0005 

IMPLICIT  REAL* 8 ( A-H , O-Z )  0006 

REAL*  4  GUN, XLOC, DATE, PRENG , ELEV , ABOVE  0007 

INTEGER  «  4  SID.S.SPOTZ  000B 

COMMON  KPOI NT , KRDN , DEFLM , DEF , RN , ALP , BET ,  BHKAX  ,  BHMI N  ,  BHBAR ,  RNG ,  TANK  0009 

1 , BH , NNTUV  0010 

DIMENSION  X ( 5  )  , Y( 5 ) , AZ ( 201 , 5 ) , ISB ( 4 ) , PT( 1 2 , 2 ) , ISC { 3 ) , CNT{ 2 )  0011 

DIMENSION  KPOINT( 201 ) , KRDN ( 201 ) , DEFLM ( 201 ) ,DEF( 201 ) ,RN( 201 )  0012 

DIMENSION  GUN ( 14 ) , XLOC ( 1 3 ) , DATE (14), PRENG (11) .ABOVE ( 6 )  0013 

DIMENSION  TIM£( 4 , 200 ) , ALPHA! 4 ,200 ) ,DIV( 200 )  0014 

DIMENSION  ALP(200) ,BET(200) ,BBH(3)  0015 

DIMENSION  BHMAX(200),BHMIN(200).BHBAR{200)  0016 

DIMENSION  RNG{ 4 ,200 ) ,TANK( 4 , 200) ,BH( 4,200 )  0017 

DIMENSION  NNTUV (  200  ) , I RND ( 5  )  0018 

DIMENSION  ALX( 5, 200  )  ,BEX( 5, 200  ) ,SID( 4  ,  200  ) , S( 50  )  0019 

DATA  SPOTZ/'S  '/  0020 

10  READ( 5 , 1040 ) IRX  0021 

GOTO  (11, 11, 12, 12, 600), IRX  0022 

12  READ( 5,1011 ) U  ,  V , HI ,HCOR  0023 

I F ( U  .E Q.  999999. VC)  GO  TO  450  0024 

121  KNTUV=KNTUV+ 1  0025 

I F ( KNTUV  .GT.  1)  GO  TO  122  0026 

DO  1200  MN-1,200  0027 

1200  NNTUV ( MN  )  =  0  0028 

122  I F ( I RX  .EO.  3)  GO  TO  17  0029 

123  I CNT« I CNT* 1  0030 

M= I CNT  0031 

L= I RND ( M  )  0032 

DO  15  K= 1 , L  0033 

I F ( ALX ( M , K ) . EO . - 9 9 9  . DO )  GO  TO  815  0034 

NNTUV ( K )« NNTUV ( K ) + 1  0035 

ALX  (  K  ,  K  )  ■  (  A.LX  (  M  ,  K  )  +  (  BEX  (  M  ,  K  )  +HCOR  )/60. DO)*. 017453293  0036 

815  B E T ( K ) « B E X ( M , K  )  0037 

ALP ( K ) « ALX ( M , K )  0038 

15  CONTINUE  0039 

17  GO  TO  ( 7 1  ,  1 0  )  ,  I  ERR  0040 

11  READ( 5, 3100  )GUN  0041 

II«1  0042 

DO  1  IQZ*=1  ,  50  0043 

1  S ( I QZ ) «  SPOTZ  0044 

READ( 5, 31 50 ) XLOC, DATE, PRENG, ELEV, ABOVE  0045 

READ(5,1000)  KRD,N,NSTEP,A  0046 

KNTUV-0  0047 

DO  13  J-1,200  0048 

BHMAX( J ) «0 . DO  0049 

BHMIN( J1-999999.D0  0050 

BHBAR( J)«0.D0  0051 

13  DIV(J)-N  0052 

IERR-1  0053 

DO  20  1-1 ,N  0054 

20  READ (5,1010)X(I),Y(I)  0055 

READ( 5,1011 )X( 5) , Y ( 5) , AZD , AZM  0056 

KDEG-AZD+. 5  0057 

DKMIN-AZM  0058 

A2 ( 1 , 5 ) « ( ( AZD+AZM/6  0.D0)-35. DO)*. 017453293  0059 

IRDG-0  0060 


1CNT-0  0061 

DO  70  J«1,N  0062 

K«1  0063 

DO  70  I«1 , 200  0064 

GO  TO  (  30 , 60  )  , K  0065 

30  READ ( 5,1030  ) A2D , AZM , ALX ( J , I ) , BEX ( 0 , 1 ) .ALPHA (0,1 ) , TIME ( J , I ) . SI D( J , I  0066 

1)  0067 

A2 (1,5) -AZ (1,5)  0068 

32  I E (AZD  .EO.  -999. DO)  GO  TO  60  0069 

33  1 P ( AZD  .NE.  999. DO)  GO  TO  50  0070 

34  IF  ( I RDG- I >35,40,40  0071 

35  IRDG=I  0072 

40  K«2  0073 

I RND (  J  ) «  I  - 1  0074 

GO  TO  60  0075 

50  A2(I,J)-( (AZD+AZM/60 .DO )-35. DO)*. 017453293  0076 

GO  TO  65  0077 

60  AZ( 1 , J )=99999.D0  0078 

65  IF(TINE(J,I  )  >70,66,70  0079 

66  DIVin-DIVl  1  )-l  .  0080 

70  CONTINUE  0081 

I RDG= I RDG- 1  0082 

LL- 1  0083 

71  DO  280  1=1 , I RDG  0084 

1F< SID< 1 , I ) . EO. S( II ) )  II-II+l  0085 

IF( S1D( 1 , I ) . EO- S( 1 I ) )  GO  TO  772  0086 

KRDNO= KRD  +  NSTEP * ( LL-1  )  0087 

LL-LL+1  0088 

772  NHR«=1  0089 

J=1  0090 

DO  90  K  = 1 , N  0091 

I F ( AZ ( I , K  )  .EO.  99999. DO)  GO  TO  90  0092 

80  NHR*NHK+ 1  0093 

1  SB ( J ) «  K  0094 

J=J+1  0095 

90  CONTINUE  0096 

GOTO  (274, 110, 120, 130, 140), NHR  0097 

110  I  SB ( 2  ) *  5  0098 

120  LSV= 1  0099 

CALL  POINT(A2,I , 1  SB , 1 ,2,X,Y,PT ,LSV)  0100 

CNT  (  1  )  *■  PT  (1,1)  0101 

CNT ( 2 ) =PT (1,2)  0102 

GO  TO  254  0103 

130  LST«4  0104 

GO  TO  150  0105 

140  LST-1  0106 

150  ARSV-1.D20  0107 

DO  240  L-LST,4  0108 

CALL  ISCB(L.ISC)  0109 

IF  (NHR-4)160,160,200  0110 

160  DO  170  1BC-1, 3  0111 

170  ISC(IBC)-ISBdBC)  0112 

200  DO  210  J-1,2  0113 

KST-J+1  0114 

DO  210  K-KST.3  0115 

210  CALL  POINT (AZ,I,ISC,J,K,X,Y,PT,L)  0116 

220  CALL  AREA( PT, L,AR )  0117 

IF  ( ARSV-AR  >240,230,230  0118 

230  ARSV-AR  0119 

LSV-L  0120 


22 


240  CONTINUE  0121 

IF  ( ARSV-A )250,250,300  0122 

250  CALL  CNTR (PT,LSV,CNT)  0123 

254  CO  TO  ( 255, 280 ) , IERR  0124 

255  GOTO  (260,270,400,400) ,IRX  0125 

260  CALL  RANGE(X(5),Y(5),AZ(1,5), CNT ,RN( I ) , DEFL , DEFLML )  0126 

KPOINT (  I  )  *  1  0127 

KRDN( I J-KRDNO  0128 

DEFLMt I )-DEFLML  0129 

DEF ( 1 ) *DEFL  0130 

GO  TO  280  0131 

270  KPOINT ( I  )  »  2  0132 

KRDN( I )-KRDNO  0133 

DEFLM ( I ) «CNT ( 1 )  0134 

DEF( I )«CNT( 2 )  0135 

GO  TO  280  0136 

274  GO  TO  ( 275, 280  )  ,  1  ERR  0137 

275  KPOl NT ( 1 ) *  3  0138 

I F ( I RX  .GT.  2)  GO  TO  280  0139 

KRDN ( I ) *KRDNO  0140 

280  CONTINUE  0141 

GO  TO  (  285 , 1  0  )  ,  IERF.  0142 

285  GOTO  (500,500,10,10,600) ,IRX  0143 

300  GOTO  (  310 , 320  ),  IERF.  0144 

310  1  ERR*  2  0145 

WRITE! 6. 4000 )  0146 

320  IF(NHR-4 >340,340,330  0147 

330  CALL  1SCR(LSV,ISC)  0148 

M«  3*LSV+1  0149 

GO  TO  350  0150 

340  M*13  0151 

350  1F(SID( 1 , I  )  ,EQ.S( II ) )  GO  TO  352  0152 

WRITE! 6 , 2031 )KRDNO,ARSV  0153 

GO  TO  353  0154 

352  IQ-II-1  0155 

WRITE(6, 2033)10, ARSV  0156 

353  DO  390  J-1,3  0157 

JISC*ISC!J)  0158 

ANUM«PT(M-J,2)-Y( JISC)  0159 

ADEN*PT(M-J,1 ) -X { JISC)  0160 

CALL  ATNSR ( ANUM , ADEN , CORK , 1 )  0161 

CORR-CORR/. 017453293+35. DO  0162 

360  I F ( CORR  .LT.  360. DO)  GO  TO  380  0163 

370  CORR* CORR- 360 . DO  0164 

380  ICORR*COKR  0165 

CORR* ( CORR-DFLOAT ( I CORR ) ) • 6 0 . DO  0166 

CORRA-CORR/60 .  0167 

CORRH- FLOAT) I CORR) +CORRA  0168 

390  WRITE (6,2032)ISC(J) ,1 CORR , CORR , CORRH  0169 

GO  TO  280  0170 

400  CALL  RANGE ( U , V, AZ (1,5), CNT , RNG( KNTUV, I ) , DEFL , DEFLML )  0171 

2F(ALP(1)  .E3.  -999. DO)  RNG( KNTUV, I )-0 .DO  0172 

IF(IRX  .EQ.  3)  GO  TO  425  0173 

KONTRL-1  0174 

IF(ALP(I)  .NE.  -999. DO)  GO  TO  420  0175 

401  TANK ( KNTUV , I ) ■ 0 .  0176 

BH ( KNTUV , I ) ■ 0 .  0177 

GO  TO  280  0178 

420  TANK (KNTUV, I )-DTAN(ALP(I ) )  0179 

BH ( KNTUV , I ) *RNG ( KNTUV , I ) *TANK ( KNTUV , I )+HI  0180 
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BHMAX (I  )  -DMAXl  (  BHMAX  (  I  )  ,  BH  (  KNTUV ,  I  )  )  0161 

BHMIN ( 1  ) -DM INI ( BHMI N ( I ) ,BH( KNTUV,  I )  )  0182 

BHBAR ( 1 ) -BHBAR ( I ) +  BH ( KNTUV , I }  0163 

GO  TO  260  0184 

425  KONTRL-2  0165 

GO  TO  280  0186 

450  I F< IERR  .EQ.  2)  GO  TO  10  0187 

DO  470  KOPY-1 , 3  0188 

WRITE! 6 , 3200 )KOPY  0189 

GO  TO  ( 460, 451 ) , KONTRL  0190 

451  DO  453  J-l, KNTUV  0191 

WRITE (6,1014)0  0192 

DO  452  I - 1 , I RDG  0193 

IF(KP01NT( I ) . EQ • 3 )  GO  TO  4530  0194 

IF(S1D( J,I ) .EQ.S(II ) )  GO  TO  4522  0195 

WRITE! 6,1016 ) KRDN ( I ) , RNG ( 0 , I )  0196 

GO  TO  452  0197 

4530  IF! SID! J , I ) .EQ. S( II ) )  WRI TE { 6 , 1 21 6 )  0198 

1 F ( SI D ( J , I ) .NE.S(II)  )  WRITE (6,1316) KRDN ( I )  0199 

1316  FORMAT ( 1 H  ,16)  0200 

1216  FORMAT ( 1 H  ,'  SPOT')  0201 

GO  TO  452  0202 

4522  WRITE (6,1116)RNG(J,I )  0203 

452  CONTINUE  0204 

453  CONTINUE  0205 

GO  TO  470  0206 

460  DO  462  I-l.IRDG  0207 

I F ( KOPY  .GT.  1)  GO  TO  461  0208 

I F ( NNTUV ( I )  .EQ.  0)  NNTUV(1)-1  0209 

BHMAX! I ) -BHMAX ( I  )-BHMIN( I )  0210 

BHBAR ( 1  ) -BHBAR ( I  ) /DFLOAT ( NNTUV  (  I  )  )  0211 

461  GOTO  (  4611, 4612, 461  3, 4614  ), KNTUV  0212 

4611  I F ( I  .EQ.  1)  WRITE! 6 ,1021 )  0213 

IF(SID( 1 ,1 ) .EQ.S! II ) )  GO  TO  4711  0214 

W^ITE! 6,1061  )  KRDN( I ) , ( ( RNG ( J , I ) .TANK! J,I ) ,BH( J, I ) ),  J- 1  ,  KNTUV ) ,  0215 

1BHBAR! I  )  .BHMAX! I  )  0216 

GO  TO  462  0217 

4711  WRITE! 6, 1161 ) ( (RNG! J, I ) .TANK! J, I ) , BH! J , I ) ), J-l  .KNTUV)  ,  BHBAR! I ) ,  0218 

lBHMAX(I)  0219 

GO  TO  462  0220 

4612  I F ( I  .EQ.  1)  WRITE! 6 , 1022 )  0221 

IF(SID( 1 ,1 ) .EQ.S( II  )  )  GO  TO  4712  0222 

WRITE! 6,1062  )  KRDN ( I ) , ( ( RNG ( J , I ) , TANK ! J , I ) , BH  (  J , I ) ) , J-l , KNTUV ) ,  0223 

lBHBARd  )  ,BHMAX(  I  )  0224 

GO  TO  462  0225 

4712  WRITE! 6, 11 62 ) < ( RNG ( J , I ) .TANK ( J , I ) , BH ( J , I ) ) . J-l .KNTUV) , BHBAR! I ) ,  0226 

1BHMAX ( I )  0227 

GO  TO  462  0226 

4613  IF! I  .EQ.  1)  WRITE! 6 , 1 02 3 )  0229  ' 

I F ( SI D ( 1 , I ) .EQ.Sfll ) )  GO  TO  4713  0230 

WRITE! 6,1063  )  KRDN!  I  )  ,  !  (  RNG ( J  ,  I ) , TANK ( J , I ) , BH! J , I ) ) , J-l .KNTUV) ,  0231 

lBHBARd  )  .BHMAX!  I )  0232 

GO  TO  462  0233 

4713  WRITE! 6,1163 )! (RNG( J, I  )  .TANK! J, I ) ,BH( J,I ) ) , J-l .KNTUV) .BHBAR! I ) ,  0234 

1 BHMAX! I )  0235 

GO  TO  462  0236 

4614  irn  .EQ.  1)  WRITE! 6,1024  )  0237 

ir(SID(l,I).EQ.S(II))  GO  TO  4714  0236 

WRITE (6, 106 4  )  KRDN! I ) , ( ( RNG ( J , 1 ) , TANK ( J , I ) ,BH ( J , I ) ) , J-l , KNTUV ) ,  0239 

1BHBAR! I ) .BHMAX! I )  0240 
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4714 

462 


A  65 

468 

469 

470 

500 


501 

503 

504 

510 

514 

515 

5999 

5151 

5551 

5152 


5352 

5252 

5452 

5153 

5353 


GO  TO  462 

WRITE ( 6 ,1164  )  (  ( RNG ( J , 1 ) , TANK ( J , I  )  ,  BH ( J , I ) ) , J-l , KNTUV) , BHBAR ( X ) , 

2  BHMAX ( I  ) 

CONTINUE 
WRITE! 6, 1017 ) 

DO  469  J-l.IRDG 

IF(SID(l,J).EO.S(II>>  GO  TO  468 
WRITE! 6 ,1019) KRDN ( J ) , BHBAR ( J ) ,BHMAX( J  ) 

GO  TO  469 

WR I TE ( 6 , 1 1 1 9 ) BHBAR ( 0 ) , BHMAX ( J ) 

CONTINUE 
CONTINUE 
GO  TO  10 

LAST-KRD+NSTEP* ( 1RDG-1 ) 

DO  560  KOPV'^l  ,  3 
WRITE (6,32CO)KOPY 

WRITE! 6 , 3300 ) GUN , XLOC , DATE , PRENG , ELEV , ABOVE 
WRITE! 6 , 3400 )KRD, LAST, KDEG.DKMIN 
IF! 1RX-2 ) 5  0  3 , 501 , 503 
WRITE!  6  ,  3700  ) 

GO  TO  504 
WRITE! 6 , 3500 ] 

KRI TE ( 6 , 36  00  ) 

DO  560  N= 1 ,1 RDG 
IBRNCH=F.  POINT! M  ) 

GO  TO  (  510, 520, 530  ,  51  0, 510  ) , IBRNCH 
Ff'EAN*  0  .  DO 
DO  515  K«  1  ,  N 
I F( DIV(M ) >515,514,515 
D I V  (  M  i  *  1  . 

PMEAN- PKEAN*TI ME! K.M)/DIV(M) 

KK-0 

DO  5999  MN-l.N 

I F ( AZ ( n , MM ) . EQ- 99999. DO)  KK-KK+l 
CONTINUE 

GO  TO  (5151, 5152, 5153, 5154), N 
IFtSID! 1 ,M) .EO. S( 1 1 ) )  GO  TO  5551 

WRITE! 6,2  001  ) KRDN ( M  )  , RN ( M ) , ALPHA ( 1 , M ) , TIME  (  1  ,  M )  ,  PMEAN 
GO  TO  560 

WRITE! 6. 2101 )RN(M) .ALPHA! 1 ,M) .TIME! 1 , M ) , PMEAN 
GO  TO  560 

I F ( KK . EO . 1 )  GO  TO  5252 

IF(SID( 1 ,M! .EQ.Stll ) )  GO  TO  5352 

WRITE(6,2002)KRDN(M)  ,  DEFLM  (  M  )  ,  DEF  ( M  )  ,  RN  (  M  )  ,AI,PHA(  1 ,  H ) ,  TIME  ( 1 ,  M  ) 

1, TIME l 2, M), PMEAN 
GO  TO  560 

WRITE (6, 2102 )DEFLM(M)  ,DEF(M) , FN{ M ) .ALPHA! 1,M) , TIME! 1 ,M ) , TIME ( 2  ,  M  )  , 
1 PMEAN 
GO  TO  560 

IF(SID(1,M).EQ.S(II) )  GO  TO  5452 

WRITE  (6,2012) KRDN ( M ) , RN! M) .ALPHA! 1 ,M) .TIME! 1 ,M ) .TIME! 2 ,H) , PMEAN 
GO  TO  560 

WRITE!  6, 2112 )RN(M>, ALPHA! 1,M) .TIME ( 1 , H ) .TIME! 2 ,M ), PMEAN 
GO  TO  560 

1F( KK . EO- 2 )  GO  TO  5253 
1F(SXD(1,M) .EQ.Stll ) )  GO  TO  5353 

WRITE  (6,2003 )KRDN(M) , DEFLM (M) , DEF! M ) , RN! M ) .ALPHA! 1 ,M) .TIME! 1 ,M ) 

1,  TIME! 2, M), TIME! 3, M), PMEAN 
GO  TO  560 

WRITE!  6, 2103 )DEFLM( M) ,DEF(M) ,RN(M) .ALPHA! 1 ,M) , TIME! l.B). TIME! 2.M), 
1TIM£(3,M) .PMEAN 
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5253 


5453 

5154 


5354 


5254 


5454 


520 


5201 

530 


5301 

560 

€00 

2110 

2114 

2014 

2102 

2002 

2100 

2000 

2010 

1000 

1010 

1011 

1015 

1016 
1116 
1017 
1119 
1019 
1014 
1021 


GO  TO  560  0302 

IF(SID( 1 ,M) .EQ.S( II ) )  GO  TO  5453  0303 

WRITE!  6 , 2  01  3  )  KRDN  (M),RN(M),ALPHA(1,M)  ,  TIME  ( 1  ,  M  )  ,  TIME  {  2  ,  M )  ,  TIME  {  3  ,  M  0304 

1). PMEAN  0305 

GO  TO  560  0306 

WRITE (6 , 2113 )RN(M) , ALPHA  (1 ,M)  ,  TI  ME ( 1 , M ) , TI ME  (  2  ,  M )  ,  TI  ME ( 3 , M ) , PMEAN  0307 

GO  TO  560  0306 

I F { K K . EQ . 3 )  GO  TO  5254  0309 

IF(SID( 1 ,M) .EQ.S(II ) )  GO  TO  5354  0310 

WRITE  (  6 , 2000  )KRDN(M)  ,  DEFLM  (  M  )  , DEF ( M !  ,  RN  (  M  )  ,  ALPHA  ( 1  ,M )  ,TIME(1 ,M)  0311 

1,TIME(2,M) ,TJME(3,M),TIME(4,M),PMEAN  0312 

GO  TO  560  0313 

WRITE!  6 ,2100)  DEFLM  (  M  )  ,  DEF  (  M  )  , RN(M)  ,  ALPHA  (1  ,  M )  .  TIME  ( 1 ,  M  )  ,  TIME  (  2  ,  M  )  ,  0314 

1 TI ME ( 3  ,  M  ) ,TIME( 4  ,  M  ) , PMEAN  0315 

GO  TO  560  0316 

I F( 5 1 0 ( 1 ,M) . EQ . S l II ) )  GO  TO  5454  0317 

WRITE!  6. 2  014 ) KRDN(M) ,RN(M) .ALPHA! 1 .M)  . TIME ( 1 , M ) , TIME! 2 , M ) .TIME ( 3 , M  0318 

1  )  .TIME!  4  ,M)  , PMEAN'  0319 

GO  TO  560  0320 

WRITE!  6 ,2114  )RN(M)  .ALPHA!  1  ,  M )  ,  TIME  ( 1  ,  M  )  .TIME  (  2  ,  M  )  .TIME  (  3  ,  M  )  ,  TIME  (  4  0321 

1 ,  M  )  ,  PJ'E.AN  0322 

GO  TO  560  0323 

IF!  SID! 1 ,K) ,EO. S( II  )  )  GO  TO  5201  0324 

WRITE!  6 , 201  0  )KKVK(K)  ,  DEFLM  (  K  }  ,  DEF(  J!)  0325 

GO  TO  560  0326 

WRITE!  6 , 211 0 ) DEFLM (K ) ,DEF(M )  0  327 

GO  TO  560  0328 

IF! SID! 1 , M ) . EQ. S! 1 1  )  )  GO  TO  5301  0329 

WRITE! 6 , 2020 ) K  R  D ” ( H )  0330 

GO  TO  560  0331 

WRITE! 6 , 2120  )  0332 

CONTINUE  0333 

GO  TO  10  0334 

CALL  EXIT  0335 

FORMAT! 1H0 ,  '  SPOT ' , 1 X , Fl 7 . 2 , IX , FI  3 . 2/ )  0336 

FORMAT!  1H  SPOT’,T14,’ - ’.T29,’ - ’  ,  T4  2  ,  F8 . 1  ,  4X ,  AB  ,  3X  ,  4  (  0337 

1F6 . 2 , 4/ )  ,  F6 . 2  )  0338 

FORMAT  ( 1 H  ,  I6.T14  ,  * - '  .T29,  ’ - '  ,  ~4  2  ,  FB  .  1 ,  4X  ,  A8 . 3X ,  4  (  F6 . 2 , 4  0339 

IX ) , F6 . 2  )  0340 

FORMAT ! 1 H  SPOT'  ,  TI  4  ,  F6 . 1 , T2 9 , F6 . 1 , T4 2 , F8 . 1 , 4X . AB , 3X , 2 ( F6 . 2 , 4X )  0341 

1.20X.F6.2)  0342 

FORMAT! 1H  .I6.T14.F6.1 .T29.F6.1 . T4 2 , F6 . 1 , 4X , A8 . 3X , 2 ( F6 . 2 , 4X ) ,20X,F  0343 

16.2)  ,  0344 

FORMAT! 1H  SPOT'  ,  TI  4  ,  F6 . 1 , T2 9 , F6 . 1 , T4 2 , FB . 1 , 4X , A8 , 3X , 5 ( F6 . 2 , 4X )  0345 

1)  0346 

FORMAT! 1H  ,I6,T14,F6.1,T29,F6.1,T42,FB.1,4X,A6,3X,5(F6.2,4X))  0347 

FORMAT! 1H0, 16 ,1X,F17.2,1X,F13.2/)  0348 

FORMAT ( 3 1 4 , F20 . 0  )  0349 

FORMAT! 2F11 .3)  0350 

FORMAT! 4F11 . 3  )  0351 

FORMAT! F4.0.1X.F3.0)  0352 

FORMAT! 1H  ,I6,F17.1,1X,F9.5,1X,F5.1)  0353 

FORMAT! 1H  SPOT* ,F17.1flXfF9.5,lX,F5.1)  0354 

FORMAT! 1H1 , IX, 'ROUND' ,9X, 'AVBHT' , 6X , ' DI  FF  '  /,  '  ’  , 1 5X , ' METERS '  ,/)  0355 

FORMAT! 1H  ,'  SPOT' , BX , F6 . 1 , 6X , F5 . 1 )  0356 

FORMAT! 1H  ,I6,8X,F6.1,6X,F5.1)  0357 

FORMAT! 1H  , 'TOWER  ',11/)  0358 

FORMAT !1H0,T24,'1'/  0359 

11H  ,'  ROUND ',1!'  RANGE  TANGENT  BHT' ) , '  AVBHT  DIFF' )  0360 
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1022  FORMAT !1H0,T24,'1',T49,'2'/  0361 

21H  ,  '  ROUND',2(’  RANGE  TANGENT  BHT'),'  AVBHT  DIFF’)  0362 

1023  FORMAT! 1H0.T24 , '1 ' ,T49, '2' ,T74 , ' 3’/  0363 

31H  ,  '  ROUND', 3('  RANGE  TANGENT  BHT’),'  AVBHT  DIFF')  0364 

1024  FORMAT (1H0,T24, '1 ’ ,T49, '2' ,  T7  4  , '3' ,T99,  '4'/  0365 

4 1 H  ROUND ' , 4 ( '  RANGE  TANGENT  BHT'),'  AVBHT  DIFF')  0366 

1030  FORMAT (2F7.2,2F7.2,A8,F6.2,A4 )  0367 

1040  FORMAT (12)  0368 


1050  FORMAT ( 1 H  16.D16.9) 

1061  FORMAT ( 1 H  ,  I  6 , 1 ( 1 X  ,  F8 . 1 , 1 X , F9 . 5 , 1 X , F6 . 1 ) , 1 X , F6 . 1 , 1 X ,  F5 . 1  ) 

1161  FORMAT ( 1 H  .'  SPOT' , 1 ( IX, F8 . 1 , 1X,F9 . 6, IX, F6 . 1 ) , IX, F6 . 1 ,1X,F5 . 1 ) 

1062  FORMAT ( 1  H  , 16,2! IX, F8.1 ,1X, F9 . 5,1X, F6 .1 ) ,1X,F6. 1 , 1X,F5.1 ) 

1162  FORMAT ( 1 H  ,'  SPOT '  ,  2 ( IX , F8 . 1 , IX , F9 . 5 , IX , F6 . 1 ) , IX , F6 . 1 , IX , F5 . 1 ) 

1063  FORMAT ( 1 H  , I  6 . 3 ( 1 X , F8 . 1 , 1 X , F9 . 5 , IX , F6 . 1 ) , 1 X , F6 . 1 , 1 X , F5 . 1 ) 

1163  FORMAT ( 1 H  ,'  SPOT ' , 3 ( IX , F B . 1 , 1 X , F 9 . S , 1 X , F6 . 1 ) , 1 X , F6 . 1 , 1 X , FS . 1 ) 

1064  FORMAT ( 1 H  , I  6 , 4 ( 1 X , F8 . 1 , 1 X , F9 . 5 , 1 X , F6 . 1 ) , 1 X , F6 . 1 , 1 X  ,  F5 . 1 ) 

2001  FORMAT  (  1 H  ,  I  6  ,  T1  4  ,  ' - ',T29,' - '  ,  T4  2  ,  F8 . 1  ,  4X  ,  A8 , 3X  ,  F6 . 2 , 34X 

1.F6.2) 

1164  FORMAT ( 1 H  ,'  SPOT  '  ,  4 ( 1 X , F8 . 1 , IX , F9 . 5 , 1 X , F6 . 1 ) , IX , F6 . 1 , IX  ,  F 5 . 1 ) 

2101  FORMAT  (  1  H  ,'  SPOT '  , Tl 4  ,  ' - ',T29,' - '  ,  T4  2  ,  F8 . 1  ,  4X  ,  A8 , 3X  ,  F6 

1.1,34X,F6.2) 

2012  FORMAT!  1H  ,I6,T14,  ' - '  ,T29,  ' - ’  ,  T4  2  ,  F8 . 1 ,  4X  ,  A8 , 3X  ,  2  (  F6 . 2 , 4 

IX) , 2 0 X , F 6 . 2 ) 

2112  FORMAT  (  1 H  ,'  SPOT '  ,  Tl  4  ,  ' - ',T29,' - '  .  T4  2  ,  F8 . 1 , 4  X  ,  A8 , 3X  ,  2  ( 

1F6.2.4X) , 2 0 X , F 6 . 2  ) 

2  013  FORMAT!  1H  ,  16  .  Tl  4  ,  ' - '  ,T2  9  ,  ' - '  ,  T4  2  ,  F8 . 1  ,  4X ,  A8 , 3X  ,  3  (  F6 . 2 , 4 

1 X  )  ,  1  0  X  ,  F  6 . 2  ) 

2113  FORMAT  (  1 H  ,'  SPOT',T14,’ - '.T29,' - '  ,  T4  2  ,  FB  .  1  ,  4X  ,  A8 , 3X  ,  3  ( 

1F6.2.4X ) , 10X.F6.2) 

2103  FORMAT ( 1 H  ,'  SPOT ' , Tl 4 , F6 . 1 , T2 9 , F6 . 1 , T4 2 , F6 . 1 , 4X , A8 , 3X , 3  (  F6 . 2 , 4X ) 
1 , 1  0  X  ,  F  6 . 2  ) 

2003  FORMAT! 1H  , I  6  ,  Tl 4  ,  F6 . 1  , T2 9 , F6 . 1 , T4 2 , F8 . 1 , 4X , A8 , 3X , 3 ( F6 . 2 , 4X )  ,  1  OX  ,  F 
16.2) 

2120  FORMAT! 1H0,'  SPOT') 

2020  FORMAT! 1H0, 16/) 

2031  FORMAT! 1H0 , 16 , IX ,D16 . 9 ) 

2032  FORMAT! 1H  ,2I4,1X,F7.3,1X,F8.3) 

2033  FORMAT ( 1H0  ,  '  SPOT ’ , 1 2 , IX , D1 6 . 9 ) 

3100  FORMAT! 14A4) 

3150  FORMAT! 1 3A4/1 4A4/1 1A4/F6 . 2 , IX , 6A4 ) 

3200  FORMAT! 1H1 , SOX, 'RANGE  FIRING  SUMMAR Y ' , Tl 1 0 , ’ COPY  ',Il,'OF  3') 

3300  FORMAT! 1  HO , ’GUN  ',14A4,'  LOCATION  ',13A4/1H  , ’ DATE  1 4A4 ,' PROJECT 
1  ENGINEER  '.11A4/1H  .’ELEVATION  OF  TRUNNIONS  ' , F6 . 2 , '  FT.  ABOVE 
26A4  ) 

3400  FORMAT! 1H0.T51, 'AZIMUTH  LINE  OF  FIR£'/1H  ,T44,'RDS.  ',16,'  TO  ',16 
1, '  AZ .  ' ,13,  '  D  ' ,F5.2, '  M  ' ) 

3500  FORMATOHO, 'RD. NO.  '  ,T18, 'DEFLECTIONS' ,T43, 'RANGE' ,T53, 'ORD  OF  DET' 
1,T85, 'TIMES  OF  FLIGHT') 

3600  FORMAT! 1H  , Tl 5 , 'MILS ', T29 ,' METERS ’ ,T4 3 , 'METERS ' ,T68 ,' 1 ' ,T78  ,' 2  '  , 
1T88, ' 3' , T9B , ' 4 ' , Tl 06 , ' MEAN ’ / ) 

3700  FORMAT! 1H0.T4 3, 'GRID  COORDINATES  OF  POINT  OF  IMPACT'/lHO, 

1'  ROUND’, 11X,'  X  '  ,  8X , '  Y  '/) 

4000  FORMAT! 1H1 , 'RANGE  FIRING  SUMMARY  ERROR  CORRECTIONS') 

END 


0369 

0370 

0371 

0372 

0373 

0374 

0375 

0376 

0377 

0378 

0379 

0380 

0381 

0382 

0383 

0384 

0385 

0386 

0387 

03B8 

0389 

0390 

0391 

0392 

0393 

0394 

0395 

0396 

0397 

0398 

0399 

0400 

0401 

0402 

0403 

0404 

0405 

0406 

0407 

0408 

0409 

0410 

0411 

0412 

0413 

0414 


SUBROUTINE  ISCR(L.l) 
IMPLICIT  REAL*8(A-H,0-Z) 
DIMENSION  1(3) 

K«  1 
DO  20 

IF  (J-L)10,20,10 
10  I  (  K  )  ■  J 

K-K+l 

20  CONTINUE 

RETURN 
END 


0415 

0416 

0417 

0418 

0419 

0420 

0421 

0422 

0423 

0424 

0425 


28 


SUBROUTINE  POI NT ( A2 , I ,15C,0,K,X,Y,PT,L)  0426 
IMPLICIT  REAL*8(A-H,0-Z)  0427 
DIMENSION  A2( 201 , 5) , ISC( 3 ) ,X{ 5) ,Y{ 5) ,PT(12 ,2 )  0428 
M” 3*L+J+K-5  0429 
UISC«1SC(J)  0430 
KISC-ISC(K)  0431 
A*=DTAN(AZ(  I  ,  JISC)  )  0432 
B«DTAN(A2( I , K I  SC ) )  0433 
PT(M,1  > - ( X ( 0 1  SC ) *A- V { JISC)  +Y(K1SC)-X(KISC)*B)/(A-B)  0434 
PT(M,2)=(PT(H,1)-X(KISC))*B+Y(KISC)  0435 
RETURN  0436 
END  0437 
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SUBROUTINE  AREA (  P  ,  L  ,  A  ) 

IMPLICIT  REAL*8(A-H,0-Z ) 

DIMENSION  P ! 1 2 , 2  ) 

I«3*L-2 
J  =  I  +  1 

A* DABS (.5D0*((P(I/2)-P(J,2))*P(K,l)+(P{U,2)-P(K,2))*P(I,l)+ 
1  (  P( K, 2 )-P( I , 2  )  )*P( J.l ) ) ) 

RETURN 

END 


0438 
0439 
0440 
0441 
04  4  2 
0443 
0444 
0445 
0446 
0447 
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SUBROUTINE  CNTR(P.L.C)  0448 

IMPLICIT  FEAL*8{A-H,0-Z)  0449 

DIMENSION  PI  12 , 2 ) ,C( 2  )  0450 

M=3*L-2  0451 

DO  1  0  1  =  1  ,  2  0452 

10  C(  I  )«  (  PI  II,  1  )+P(  M+l  ,  I  )  +P  (  M  +  2 , 1  )  )/3  .DO  04  53  * 

RETURN  0454 

END  0455 
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SUBROUTINE  RANGE ( X , Y , AZ , CNT , RNG , DEFL , DEFLM )  0456 

IMPLICIT  REAL*8{A-H,0-Z )  0457 

DIMENSION  CNT ( 2 )  0458 

RNG=  CNT (  1  ) - X  0459 

DEFL«CNT( 2 )-Y  0460 

CALL  ATNSR ( DE FL , RNG , BETA , 1 )  0461 

DEFLM=BETA-AZ  0462 

10  IF  <DEFLM-*  3. 141  59265  )  30, 30, 20  0463 

20  IF  (DEFLM-3  .  1  41  5926  5  )  50, 40, 40  0464 

30  DEFLM=DEFLM+6 . 2831 8531  0465 

GO  TO  10  0466 

40  DEFLM-DEFLM-6. 28318531  0467 

GO  TO  10  0468 

50  RNG-RNG/DCOS ( BETA )  0469 

DEFL«RNG*DSIN( DEFLM)  0470 

DE  FLM  = 1018.59164*  DE  FLM  0471 

RETURN  0472 

END  0473 
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SUBROUTINE  ATNSR ( A , B , C , L ) 

0474 

IMPLICIT  REAL*8(A-H,0-Z) 

0475 

P<=  3 . 1  4  1  5926  5 

0476 

C-DATANf A/B ) 

0477 

IF  (  A  )  1  0 , 4  0 , 4  0 

0478 

10 

IF  (B)20, 30,30 

0479 

20 

GO  TO  ( 21 , 22 ) , L 

0480 

21 

OP+C 

0481 

GO  TO  50 

0482 

22 

C«C-2 .DO’F 

0483 

GO  TO  50 

0484 

30 

GO  TO  (  31  , 5  0 ) , L 

0485 

31 

C«2.D0«P+C 

0486 

GO  TO  50 

0487 

40 

IF  (B  )20, 50, 50 

0488 

50 

RETURN 

0489 

END 

0490 
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APPENDIX  C: 

SAMPLE  INPUT  TO  TRIANGULATION  COMPUTER  PROGRAM 
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INTENTIONALLY  LEFT  BLANK. 
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01 

1  05mm 
BALL  PG. 

22  MAY  85 
SCHNELL 
13.89  MLW 


11 

4 

1  50. 

3657 

.607 

12601.626 

3868 

.530 

13870.487 

3754 

.335 

13094 . 383 

2883 

.79 

11990 . 31 

3812 

.273 

12475.603 

34. 

22  . 

18  . 

34  . 

23  . 

17. 

34. 

25. 

18  . 

34  . 

18  . 

17. 

34. 

11  . 

5. 

34. 

17  . 

5  . 

34  . 

05. 

5. 

3  4  . 

17  . 

5. 

•+999 

11 . 

08. 

1  9  . 

o  _ 

5  . 

1  6  . 

11. 

1  3  . 

19. 

9. 

42  . 

17  . 

1C. 

52  . 

6  . 

11. 

19  . 

6  . 

10. 

36. 

6  . 

11  . 

03  . 

5 . 

+  999. 

26. 

57  . 

20. 

27. 

28  . 

18. 

27. 

52  . 

19. 

27. 

24  . 

16. 

27. 

42  . 

6  . 

27. 

57  . 

5. 

27. 

32  . 

6  . 

27. 

45. 

5. 

+  999  . 

47. 

37  . 

14  . 

47. 

55. 

13. 

47. 

35. 

14  . 

4e. 

03. 

13. 

47. 

26  . 

4  . 

47. 

24  . 

4  . 

47. 

33. 

4  . 

47. 

35. 

4  . 

+  999. 

04 

3657. 

607 

12801.626 

04 

3868. 

530 

13870.487 

04 

3754. 

335 

13094 . 383 

04 

2883. 

79 

11998.31 

04 

999999. 

05 

41.  26. 


51  . 

HOAB 

15.67 

45. 

ILLUM 

16.11 

46  . 

HOAB 

15.77 

41  . 

I  LLUM 

15.94 

35. 

HOAB 

12.28 

09. 

GREN 

12.02 

3  5 . 

GREN 

11.46 

08  . 

GREN 

11.72 

07  . 

16.03 

01  . 

16.39 

05  . 

15.80 

42  . 

09  . 

11.55 

37  . 

11.63 

01  . 

11  .69 

32. 

11.99 

04  . 

50. 

45. 

40. 

45. 

55. 

35. 

50. 

48. 

15.70 

38. 

17.40 

56. 

15.73 

41  . 

17.90 

35. 

11.52 

19. 

12.07 

39. 

11.40 

17. 

11.76 

32.13 
13.3 

4.77 

12.13 
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APPENDIX  D: 

SAMPLE  OUTPUT  FROM  TRIANGULATION  COMPUTER  PROGRAM 
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Intentionally  left  blank. 
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^  o- 
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ph 

rw 

pH 

pH 
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ph 

•H 

•u 

pH 

pH 

pH 
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s 
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